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Figure  1:  A  feedback  model  for  the  Colpitts  or  Hartley  oscillator  used  to  drive  rf-discharge  lamps  in  vapor-cell 
atomic  clocks.  For  the  Colpitts  oscillator,  Z2  and  Z3  are  capacitors  with  Zj  an  inductor.  For  the  Hartley  oscillator, 
Z2  and  Z3  are  inductors  with  Z\  a  capacitor.  In  either  case,  an  inductor  is  assumed  to  surround  the  glass  bulb  of 
the  rf-discharge  lamp,  thereby  providing  energy  to  ionize  the  Rb  atoms  and  generate  resonant  light. 


Abstract 

In  the  rf-discharge  lamp  of  an  atomic  clock ,  the 
inductor  of  a  Colpitts  or  Hartley  oscillator  surrounds  a 
glass  bulb  containing  a  vapor  of  Rb  and  a  noble  gas 
(typically  Xe  or  Kr).  Rf-energy  is  extracted  from  the 
field  leading  to  ionization  of  the  Rb,  and  in 
recombination  with  electrons  these  Rb  ions  produce  the 
resonant  light  necessary  for  atomic  signal  generation. 
From  an  electrical  perspective,  the  discharge  can  be 
viewed  as  a  permeable  medium  located  inside  an 
inductor's  coils .  This  permeable  medium,  however, 
must  have  both  a  real  and  an  imaginary  part:  not  only 
does  the  discharge  alter  the  phase  of  the  circuit's  rf- 
field,  it  also  extracts  energy  from  the  resonant  circuit. 
Here,  we  consider  the  manner  in  which  this  complex 
permeability  enters  the  electrical  description  of  the 
oscillator,  and  its  likely  dependence  on  discharge 
parameters. 


I.  Introduction 

We  consider  the  very  general  oscillator  feedback 
circuit  illustrated  in  Fig.  1,  where  an  external  AC 
voltage  V'i  is  added  to  a  feedback  signal,  Vf,  to  produce 
the  input  voltage  Vj  for  an  amplifier.  Of  course,  the  key 
element  in  the  figure  is  the  network  of  complex 
impedances  in  the  feedback  loop,  which  typically  take 
one  of  two  configurations  [1]:  a  Colpitts  oscillator 
configuration,  where  Z2  and  Z3  are  capacitors  with  Zi 


the  loop-inductor  surrounding  a  lamp’s  glass  bulb;  or  a 
Hartley  oscillator,  where  Zx  is  a  capacitor,  and  Z2  and 
Z3  are  inductors  (one  of  which  corresponds  to  the  wire 
loops  surrounding  a  lamps’  glass  bulb  [2,3]). 


Figure  2:  Figure  1  redrawn.  Here,  Vf  is  the  input 
voltage  to  the  amplifier. 

Focusing  for  the  moment  on  the  feedback  portion 
of  Fig.  1,  this  can  be  redrawn  as  Fig.  2,  and  from  this 
figure  it  is  straightforward  to  show  that 
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V, 


vfz,z2 


0  -  r,(z1  +  z2  +  z3)+zi(z2+z3)  ~  K 

Then,  returning  to  Fig.  1 ,  we  see  that 

V0  =gvi  =  g(v;+vf), 

and  from  Eq.  (1)  this  yields 

gv; 


=  — L*  (1) 


V  =- 


1-gK 


(2a) 


(2b) 


The  only  way  V0  can  be  non-zero  without  an  input 
signal  V’,  (i.e.,  the  only  way  the  circuit  can  self- 
oscillate)  is  if  gK  =  1.  Since  g  is  real,  self-oscillation 
implies  that  k  must  also  be  real.  In  other  words,  the 
feedback  signal  must  return  to  the  amplifier  input  in 
phase  (i.e.,  with  Vf  =  |Vfje10,  0  must  equal  2n7t  where  n 
is  an  integer).  Thus,  for  self-oscillation  we  require  that 

Re[x]  =  -  and  lm[tc]=0.  (3) 

g 

If,  to  first  order,  we  assume  that  the  Z,  are  pure 
inductors  and  capacitors,  then  from  Eq.  (1)  it  is  clear 
that 


Im[K]  =  Z,  +  Z2  +  Z3  =  0  •  (4a) 


this  as  a  resonant  RLC  circuit  as  illustrated  in  Fig.  3. 
Briefly,  Cs  and  Rs  are  the  series  capacitance  and 
resistance,  respectively,  which  must  exist  for  real  loops 
of  wire;  while  the  parallel  capacitor,  Cp,  represents  the 
capacitance  that  must  exist  between  the  loops  of  wire 
and  the  lamp's  metal  housing. 


Defining  Zl  as  the  complex  impedance  of  the  RLC 
circuit,  it  is  straightforward  to  show  that  in  the  absence 
of  a  discharge  (i.e.,  the  inductor  loops  simply  surround 
air)  that 
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In  these  expressions, 


co.  = 


C  +Cn  j  c  +  c 
s  p  and  f  =  s  p 


LCsCp 


(5b) 
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Re[K]  =  hlh.  ^  =  g.  (4b) 

^2  ^1 

Thus,  in  order  for  g  to  be  positive,  we  see  from  Eq. 
(4b)  that  oscillation  limits  the  choices  for  the  Z,:  if  Zj  is 
an  inductor  then  Z2  must  be  a  capacitor;  alternatively,  if 
Zj  is  a  capacitor  than  Z2  must  be  an  inductor. 

II.  The  Inductor  in  a  Real  Discharge  Lamp 

L  Cs  Rs 


cP 


Figure  3:  In  a  real  rf-discharge  lamp,  the  impedance 
Zj  that  corresponds  to  the  loops  of  wire  surrounding 
the  lamp’s  glass  bulb  should  be  replaced  with  the 
RLC  circuit  shown  here.  For  the  Colpitts  lamp 
oscillator  that  we  used  in  our  experiments  (i.e.,  L  = 

650  nH),  one  of  us  (FW)  measured  Cs  =  3.7  pF.  R,,  = 

8  Q,  and  Cp  =  9  pF. 

Recent  measurements  made  by  one  of  us  (FW) 
clearly  demonstrate  that  in  a  real  rf-discharge  lamp 
oscillator  there  is  a  subtlety  that  must  be  included  in 
the  model  of  Section  I.  Briefly,  whatever  element  in  the 
feedback  network  corresponds  to  the  loop  inductor 
surrounding  the  lamp’s  glass  bulb,  we  must  consider 


Figure  4  shows  plots  of  Re[ZL]  and  Im[ZL]  as  functions 
of  frequency  for  the  parameters  that  were  measured  for 
our  experimental  setup,  and  which  are  given  in  the 
caption  of  Fig.  3. 


Frequency,  MHz 

Figure  4:  Real  (blue)  and  imaginary  (red)  parts  of  the 
complex  impedance  for  the  RLC  model  of  the 
inductor  loops  that  surround  the  lamp's  glass  bulb 
(assuming  no  discharge).  We  employed  the 
parameters  given  in  the  caption  of  Fig.  3. 

As  illustrated  in  Fig.  5,  to  include  the  influence  of 
the  discharge  in  Eqs.  (5),  we  let  L  -»  pL,  where  p  is  a 
complex  scalar  permeability:  p  =  a  -  i(3  [4].  Here,  a 
describes  that  part  of  the  plasma  discharge  that  gives 


2 


rise  to  a  phase  shift  of  the  rf-field,  and  which  we  expect 
will  influence  the  resonant  frequency  of  the  RLC 
circuit,  coL.  p  describes  rf-energy  extraction  by  the 
discharge,  and  so  we  expect  it  to  contribute  to  Rg. 


Figure  5:  In  the  rf-discharge  lamp  of  an  atomic 
clock,  the  inductor  coils  of  the  Colpitts  or  Hartley 
oscillator  surround  a  glass  bulb  that  contains  Rb  and 
a  noble  gas  (e.g.,  Xe  or  Kr).  Since  rf-energy  is 
extracted  from  the  field  to  ionize  the  Rb,  the 
permeability  of  the  material  inside  the  inductor  must 
have  an  imaginary  part  that  leads  to  resistance. 

Without  going  into  all  the  details  here,  an  analysis 
of  Fig.  3  including  a  term  for  the  discharge’s 
permeability  gives  credence  to  these  expectations. 
Thus,  in  order  to  include  the  discharge  in  Eqs.  (5)  and 
(6)  we  need  only  make  the  replacements 


Figure  7:  Real  (blue)  and  imaginary  (red)  parts  of  the 
complex  impedance  for  the  RFC  model  of  Fig.  3  with  a  =  1.2 
and  p/a  =  0.03. 


As  illustrated  in  Figs.  5,  6,  and  7,  Im[ZL]  crosses 
zero  at  two  locations  near  coL.  From  Eq.  (5b),  it  is 
straightforward  to  calculate  these  zero  crossings  in  the 
case  that  coCsRs  «  1,  which  will  typically  be  the  case 
(even  when  the  discharge  is  present).  We  define  cos  as 
the  series  resonant  frequency  and  cop  as  the  parallel 
resonant  frequency: 


and 


=  /Cs+Cp 

L  trLC°cp 


R 


R.  +|  — 
a 


(7a) 


(7b) 


Figure  6  shows  plots  of  Re[ZL]  and  Im[ZL]  with  a 
discharge  present:  a  =  1.2  and  p/a  =  0.01.  Figure  7  is 
similar,  but  with  p/a  =  0.03. 
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Figure  6:  Real  (blue)  and  imaginary  (red)  parts  of  the 
complex  impedance  for  the  RLC  model  of  Fig.  3  w  ith  a  =  1.2 
and  p/a  =  0.01. 
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Here,  the  plus  and  minus  signs  refer  to  cos  and  cop, 
respectively,  and  we  note  from  Eq.  (6)  that  f  >  1. 
Typically,  the  oscillator  will  resonate  at  cos,  since  it  is 
at  this  frequency  that  the  real  part  of  the  impedance  is 
minimized. 


III.  The  Colpitts  Oscillator 

As  a  consequence  of  the  considerations  presented  in 
Section  II,  and  referring  to  Fig.  2,  for  a  Colpitts 
oscillator  we  take 


Z3  —  Rc 


coC, 


Z,  =  Re[ZL]+iIm[zJ. 


(9a) 

(9b) 


Here,  Rc  is  the  equivalent  series  resistance  that  we 
expect  for  a  real  capacitor,  but  which  we  will 
henceforth  assume  is  zero. 


To  determine  the  resonant  frequency  of  the  Colpitts 
oscillator  we  employ  Eq.  (4a)  along  with  Eqs.  (5),  (7) 
and  (8).  First,  however,  we  note  that  in  Colpitts 
oscillator  circuit  used  in  our  experiments  we  have  C3  = 
150  pF  and  C2  =  30  pF.  Thus,  if  we  define  a  “Colpitts 
capacitance,”  CCP,  as 
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c 


CP  ~ 


c2c3 
C2  +C3 


(10) 


then  in  our  case  we  have  Cc p  =  25  pF,  so  that  Cp/Ccp  = 
0.36  and  Cs/Ccp  =  0.15.  Thus,  while  the  Colpitts 
capacitance  is  greater  than  the  capacitances  occurring 
in  the  RLC  circuit  of  the  lamp’s  coils,  it  is  not 
significantly  greater.  Consequently,  the  complex 
impedance  of  these  capacitors  cannot  be  ignored,  and 
we  have  from  Eq.  (4a) 
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Figure  8:  The  parallel  resonant  frequency  (red),  cop, 
and  series  resonant  frequency  (blue),  os,  of  the 
Colpitts  oscillator. 
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Figure  9:  The  difference  between  the  parallel 
resonant  frequency  (red)  and  the  series  resonant 
frequency  (blue)  from  their  values  in  the  case  of  Ccp 
— »  oo.  Notice  that  the  units  for  the  parallel  resonant 
frequency  difference  are  kilohertz. 


Solving  Eq.  (11)  numerically  for  the  series  and 
parallel  resonant  frequencies,  cos  and  cop,  respectively. 
Figs.  8  and  9  show  cos  and  cop  as  functions  of  Ccp  with 
a  =  1.2  and  p/a  =  0.01.  Figure  8  shows  the  actual 
resonant  frequencies,  while  Fig.  9  shows  the  difference 
between  these  frequencies  and  the  resonant  frequencies 
in  the  case  that  CCp  ->  oo.  (Notice  in  Fig.  9  that  the 
units  for  Acop  are  kHz.)  Clearly,  the  Colpitts 
capacitance  has  a  small  but  non-negligible  effect  on  cos, 
affecting  its  value  by  -  4  %  in  our  case  of  Ccp  =  25  pF. 
Alternatively,  the  Colpitts  capacitance  has  hardly  any 
effect  on  cop,  affecting  its  value  by  -  0.016  %  in  our 
case.  Thus,  while  the  Colpitts  capacitors  C2  and  C3 
cannot  be  ignored,  the  resonant  frequency  of  the 
oscillator  is  primarily  determined  by  the  RLC  circuit  of 
the  lamp  coils  illustrated  in  Fig.  3.  This  will  most 
certainly  vary  from  lamp  to  lamp  (even  for  the  same 
circuit  design),  and  likely  gives  rise  to  some  of  the 
variability  among  Rb  clock  lamp  oscillators. 

IV.  Physical  Model  for  p 

In  this  section,  we  attempt  to  tie  p  to  physical 
characteristics  of  the  discharge.  To  begin,  we  first  note 
that  as  a  resistance  term,  we  expect  P  to  be  proportional 
to  the  circuit’s  power  loss.  More  specifically,  we 
expect  p  to  represent  electrical  power  that  flows  out  of 
the  circuit,  into  the  discharge,  and  from  the  discharge  is 
irreversibly  lost . 

Without  too  much  difficulty,  we  can  imagine  at 
least  three  processes  that  lead  to  irreversible  energy 
flow  out  of  the  discharge:  discharge  heating  of  the 
bulb’s  glass  walls  [5,6],  electron  excitation  of  Xe  and 
the  resulting  Xe  photon  emission  [7],  and  electron/Rb+ 
recombination  leading  to  photon  emission  [7].  We  do 
not  include  Rb  ionization  in  this  list,  since  ionization 
by  itself  does  not  represent  energy  loss;  energy  is  only 
lost  by  the  discharge  when  those  ions  recombine  with 
electrons  and  emit  a  photon  that  escapes  the  discharge. 
In  this  regard,  it  is  important  to  note  that  radiation 
trapping  [8,9]  likely  limits  the  discharge’s  energy  loss 
by  electron/Rb+  recombination,  since  the  energy 
carried  by  the  radiation-trapped  photon  has  a  high 
probability  of  getting  back  into  the  discharge. 

We  therefore  write 

P  ~  7th (t “ Tdc)+  ^X^,2  e  dx 

+  (yDn2-rRT([RbD).  (12) 

The  first  term  on  the  right-hand-side  of  Eq.  (12) 
represents  rf-heating:  TDC  is  the  temperature  of  the 
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discharge  when  the  rf-field  supplies  no  additional  heat 
(i.e.,  it  is  the  temperature  of  the  discharge  as  defined  by 
a  DC  heater  around  the  lamp  bulb)  and  1/Yth  is  a 
thermal  time  constant.  The  second  term  on  the  right- 
hand-side  of  Eq.  (12)  corresponds  to  electron  excitation 
of  Xe:  Si  is  the  first  excited  state  of  Xe  at  8.4  eV,  and 
Te  is  the  electron  temperature.  From  our  lamp’s  spectra 
we  estimate  Te  ~  3500  K  [5],  so  that  kTe  -  0.3  eV. 
Finally,  the  last  term  in  brackets  on  the  right-hand-side 
of  Eq.  (12)  represents  electron/Rb+  recombination. 
Since  we  expect  charge  neutrality  in  the  discharge,  the 
density’  of  Rb+  should  equal  the  density  of  free 
electrons,  and  the  rate  of  recombination  should  be 
proportional  to  those  two  densities  multiplied  by  an 
ambipolar  diffusion  time  constant,  l/yD.  Further,  we 
reduce  the  loss  of  electron/Rb+  recombination  by  a 
radiation-trapping  term,  T RT,  which  we  expect  will  be 
some  complicated  function  of  the  neutral  rubidium 
atom  density  in  the  discharge,  [Rb] . 


Focusing  on  the  Xe  excitation  term,  we  write 


Pxe  - 


(kTj 

Then,  evaluating  the  integral  we  get 


dx  • 


e'*,/lT*  +  — Erfc 

£1 

||  kTe  2 

JkTc. 

(13) 


(14) 


where  Erfc[...]  is  the  complimentary  error  function 
[10].  Note,  however,  that  Si/kTe  ~  28,  so  that  the  first 
term  on  the  right-hand-side  of  Eq.  (14)  has  a  value  of  ~ 
4x1  O'12,  while  the  second  term  on  the  right-hand-side 
of  Eq.  (14)  has  a  value  of  ~  6x1  O'14.  Thus,  to  good 
approximation  for  our  limited  range  of  Tc  values  we 
can  ignore  the  second  term  on  the  right-hand-side  of 
Eq.  (14)  and  write 

Pxc  =  Yxc 


Figure  10  shows  Ap/Apmax  and  ApXe/Apxe,max  as 
functions  of  lamp  temperature  from  our  measurements 
of  a  Colpitts-oscillator  lamp.  Briefly,  we  measure  the 
complex  impedance  of  the  lamp  coils  illustrated  in  Fig. 
3  at  the  oscillation  frequency,  and  from  that 
measurement  determine  a  and  p.  Here,  we  define 
Ap/Apmax  as 

AP  =  P~P-m  ,  (16) 

APmax  Pmax  P  min 

with  a  similar  expression  for  Apxe/Apxc,max-  The  fact 
that  Ap/Ap^  and  ApXe/Apxe,max  track  each  other  so 
well  suggests  that  our  physical  model  for  p  has  value. 


Additionally,  it  suggests  that  at  low  lamp  temperatures 
(i.e.,  T  <  138  °C)  changes  in  p  derive  principally  from 
changes  in  electron  temperature.  Though  we  don’t  as 
yet  have  a  good  understanding  of  why  p  increases  at 
high  temperatures  (i.e.,  T  >  138  °C),  it  may  be  that  the 
electron  density  in  the  discharge  increases  at  these 
higher  temperatures,  or  that  radiation-trapping  is  less 
efficient  at  returning  the  energy  of  electron/Rb+ 
recombination  photons  to  the  discharge  (and  hence  the 
electrical  circuit). 
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Figure  10-  Ap/Apmax  (yellow  squares)  and 
Apxe/Apxe,max  (orange  circles)  as  a  function  of  lamp 
temperature  for  our  lamp  operating  at  nominal  rf- 
power. 


V.  Summary 

In  analyzing  an  atomic  clock’s  rf-discharge  lamp 
oscillator,  we  found  that  the  lamp  coils  must  be 
modeled  as  an  RLC  circuit  if  the  electrical 
characteristics  of  the  oscillator  are  to  be  properly 
understood.  Further,  from  an  electrical  perspective,  the 
plasma  inside  the  lamp  cannot  be  ignored,  and  that  this 
should  be  included  in  the  circuit  analysis  as  a 
permeable  medium  with  a  real  and  an  imaginary  part:  p 
=  a  -  ip.  Re[p]  plays  a  primary  role  in  determining  the 
lamp  circuit’s  resonant  frequency,  while  Im[p]  plays  an 
important  role  in  determining  the  circuit’s  energy  loss 
(i.e.,  the  circuit’s  Q).  Finally,  we  considered  the 
discharge  processes  that  likely  contribute  to  Im[p],  and 
presented  evidence  for  a  primary  played  by  electron 
temperature. 
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